Emergence of Clusters in Growing Networks with 
Aging 

Nuno Crokidakis and Marcio Argollo de Menezes 

Instituto de Ffsica, Universidade Federal Fluminense 

Av. Litoranea s/n 

24210-340, Niteroi - RJ, Brazil 

E-mail: nuno@if.uff.br, marcio@if.uff.br 

Abstract. We study numerically a model of nonequilibrium networks where nodes 
and links are added at each time step with aging of nodes and connectivity- and 
age-dependent attachment of links. By varying the effects of age in the attachment 
probability we find, with numerical simulations and scaling arguments, that a giant 
cluster emerges at a first-order critical point and that the problem is in the universality 
class of one dimensional percolation. This transition is followed by a change in the 
giant cluster's topology from tree-like to quasi-linear, as inferred from measurements 
of the average shortest-path length, which scales logarithmically with system size in 
one phase and linearly in the other. 



PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 64.60.-i 



Emergence of Clusters in Growing Networks with Aging 



2 



1. INTRODUCTION 

The understanding of natural, technological and social phenomena through the network 
perspective has motivated a large body of research in the past 10 years [HOE], as many 
of their properties can be inferred from the simple idea that relations among elements 
in such systems can be interpreted, at a certain level of abstraction, as nodes and links 
of a complex network. Examples range from physically interacting proteins in the cell 
[U [5] and the set of routers comprising the Internet [B] to social networks pO, [8] and, 
according to empirical measurements of a few network properties, like degree distribution 
and correlations, clustering and average shortest path, it is observed that many of these 
networks are strikingly similar (in the statistical sense) [2]. Thus, generic features of 
these networks can be predicted by the analysis of simple models. Since the classic work 
of Erdos and Renyi 0, HD], a whole set of interesting results in percolation theory have 
been brought up revealing the importance of network structure in determining critical 
properties [TH H21 [131 El EEl [16] . In particular, it has been found that in nonequilibrium 
(growing) networks with exponential [11] or power-law degree distribution [12] an infinite 
order critical point separates a phase with many finite clusters from another where a 
single macroscopic connected cluster emerges. It has been shown both numerically [H] 
and analytically [12] that all the derivatives of the average size of the largest cluster, 
taken as the order parameter, are zero at the critical point. Here we study a similar 
problem, where networks are grown by the addition of nodes and links that attach 
preferentially, but with an age-dependent probability. With computer simulations and 
scaling arguments we show that one has a first-order transition in the size of the giant 
connected cluster as one makes it less likely to attach links to older nodes. Following the 
transition the topology of the largest cluster changes from tree-like to one-dimensional, 
as one can infer from measurements of the average shortest-path length inside the largest 
cluster in each phase. 

2. THE MODEL 

Let a network grow from an initial cluster of mo fully-connected nodes (we use tuq = 2) 
by the addition of a node and a link at each time step. This new link will randomly join 
a pair of nodes with probability H(ki,kj,ai,aj,t) = U(ki,ai,t)U(kj,aj,t), where 
Il(fc, a, t) = C(A + k)e~ a ^~ a ^ is the probability that a node with k links and added at 
time a < t to the network receives a link and C(t) = J2i^(h,ca,t). Self-links and 
multiple links between the same pair of nodes are forbidden. The parameter a plays 
the role of the inverse of a timescale /3, which suppresses the attractiveness of older 
nodes for new connections. The case a = (j3 — > oo) has been extensively studied in 
[12] , who reported an infinite-order phase transition as the number b of links added per 
time step is varied around a critical value b c . Here we analyze the effect of varying the 

| an initial attractiveness [17] Aq > is necessary for the newly added node, which has k = 0, to 
participate in the dynamics. We set Aq = 1 in the following. 
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memory parameter a and find another nontrivial phase transition at a critical value 
a c separating a phase with multiple clusters from another where a single, O(N), giant 
component emerges. 

For each value of a we grow 10 3 networks with up to iV = 15000 nodes. Since 
we do not force the newly added node to immediately attach to a preexisting cluster, 
there is always a possibility to create networks with isolated clusters. For a given value 
of a, we identify these clusters in each network with breadth- first-search and calculate 
the average size of the largest cluster (S(a)) in each network, normalized by the total 
number of nodes N (Fig. [1]). 




Figure 1. The relative size of the largest cluster (left) and the average number of 
clusters (right) as a function of a. We observe a phase transition at a given value a c 
where one giant cluster emerges. Results are averaged over 10 3 samples and clusters 
are identified with breadth- first-search. 



For small values of a the network is fragmented in many small clusters, whereas 
after a (size- dependent) value a c (N) a giant component most likely exists. In order 
to determine the value of a c (N) we analyze the fluctuations of the order parameter 
x(a, N) = (S 2 ) — (S) 2 as a function of a and identify the position of the maximum of 
X with a c (N) (Fig. [2] on the left). Finite-size scaling analysis of usual (first or second 
order) phase transitions suggests a power-law scaling for the critical point shift 

a c (N) = a c (oo) + A iV" 1 ^, (1) 

where A is a constant and v = 1/d at a first-order (discontinuous) phase-transition [19|. 
Nevertheless, we find that for this model the position of the maxima of the susceptibility- 
like parameter x(a c (iV), N) scale as 

\na c (N) =\na c (oo) - K N- 1 , (2) 

where lna c (oo) = 1.717(1) and K = 516.1(7). In the thermodynamic limit (N — > oo) 
we obtain a c = 5.568(2) (We refer to ct c (oo) as a c in the following). 

As another unusual feature of this phase transition, we find that the relative size 
of the largest cluster {S)/N is a function of the ratio a/a c (N), as can be seen in Fig. 
[3j This might result from the fact that, close to a c , the dimensionless scaling variable 
for the average size of the largest cluster (S) is x = log(a), as appears from the scaling 
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Figure 2. Fluctuation of the order parameter for different network sizes as a function 
of a (left) and the position of their maxima a c (N) versus the inverse of network size 
./V -1 (right). The straight line corresponds to the best fit In y = 1.7171 — 516.17 x. 



of the rounding of the phase transition [Eq. (El)], so a function of distance in variables 
x translates into a function of the ratio in variables a. It is noteworthy to mention 
that the same occurs in Id percolation |18j . where the characteristic length scales 
as log^/pe)^ 1 which can be expanded as £(p) ~ \p — p c | _1 when p is close enough to 
Pc = 1- 




(a/a )exp(b N ) 



Figure 3. Scaling plot of the order parameter < S > /N as a function of ■ 



7W) 



where 



a c (N) 



,-6iV" 



. The values of the parameters are a c — 5.5682 and b = 516.17. 



The analogy with percolation on a ring extends further when we look at the cluster 
size distribution (Fig. HJ). In the limit a — > 0, depending on the number b of bonds 
added per time step, there is a critical phase without a giant connected component or 
a normal phase with a giant cluster (as predicted by Mendes et al. in [12])- Since 6=1 
in our case, there is no giant cluster at a — > 0, and the network consists of clusters of 
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Figure 4. The relative cluster size distribution n(s) for some values of the parameter 
a and N = 15000 nodes. The exponential decay of n(s) for every value of a suggests 
that t = 0, as in 1— dimensional percolation. 

many different sizes s, distributed in a pure exponential form 

n(s) ~ e~ s/ ^ a) (3) 

which suggests that r = [§|. Moreover, as seen in the inset of Fig. HI the whole 
a > a c phase is critical, in the sense that £(a > a c ) = oo. One might guess this 
result by noting that when a 3> a c , the most typical configuration is a line of nodes, 
that is, each added node gets connected to its immediate neighbor in the past. Given 
P(a, b) = C(k a + l)(kb + l)e~ a ( ta+th > , the probability of joining nodes a and b with ages 
t a and tb and degrees k a and respectively, and neglecting, as a first approximation, 
terms with contributions smaller than e~ 4a , the probability that a newly added node 
gets connected to this line is P con ~ 4/27(2e~ Q + 3e~ 2a + 3e~ 3a + 3e _4a ) and a new 
cluster emerges with probability P n ~ (24/27)e _4Q . These probabilities get comparable 
when 4 < a < 5 but, since there is always a nonzero chance of not connecting the new 
node to the previous ones one finds that for large values of a the system should have 
scaling properties of a critical one-dimensional percolation network. 

We also find power-law scaling for the divergence of the susceptibility at the critical 
point 

X (a c ) ~ N\ (4) 

with 7 = 1.0(2), as depicted in Fig. [U Sufficiently close to a c , one can approximate 
the equation for the phase transition shift, Eq. (j2J), by ln(a c (oo)/a c (iV)) ~ (a c (N) — 
a c (oo))/a c (oo) ~ iV _1 , and we find that v = 1. This is what one expects at first-order 
phase transitions in one-dimensional systems, based on scaling and renormalization 
group arguments [19]. 

§ The general scaling function for the distribution of cluster sizes writes n(s) ~ s~ T e -s /£( Q ) with a 
characteristic cluster size diverging near a c as £(a — > a c ) ~ \p — p c \~ v ■ In Id one has i/ = 1 and r = 0. 
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Figure 5. Maximum of the susceptibility x(a c (iV) , N) versus N. The best fit indicates 
a power-law divergence x( a c) ~ ^ 7 with 7 = 1.0(2). 




N N 

Figure 6. Scaling of the average shortest-path between every pair of nodes in the 
largest connected component with system size N, for values of a below (left panel) 
and above the critical point a c (right panel). The emergence of a giant connected 
component is followed by a change from tree-like to quasi-one dimensional topology, 
as indicates the crossover from logarithmic to linear scaling of (I) with N. 

In the limit a — > 0, when the age-dependence becomes vanishingly small, links 
between nodes are added in a random fashion (although with different rates on nodes 
with different degrees) and the finite clusters are locally trees in the large size limit 
[T2| [15] . In this limit one should expect a logarithmic dependence of the average shortest- 
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path length (l(N,a)) with system size N 
1 N 

l = jrE^)-OHN), (5) 

where Uj is the minimum number of links that must be traversed to join nodes % and j, 
and ( ) means average over realizations. On the other side, on rf-dimensional networks 
I is proportional to the linear dimension L. Thus, as one increases the effects of aging 
in the preferential attachment of links, one expects a transition from "small" to "large- 
world" networks [20], or a change from logarithmic to linear scaling of the average 
shortest-path with system size. In Fig. [6] we show the average shortest-path (l(N,a)) 
for every pair of nodes in the largest cluster of each network generated for different 
system sizes N. For a < a c we find logarithmic scaling of (/) with N, while (/) ~ N 
for a > a c , supporting our view of a first-order transition in the universality class of Id 
percolation for this problem. 

3. Conclusions 

We have studied the percolating properties of growing networks with age and degree 
preferential attachment of links: nodes introduced earlier in time are exponentially less 
likely to acquire new links than "younger" ones and links attach preferentially to nodes 
with high degree. One node and one link are introduced at each time step and the effect 
of aging on a node with age a and k links is varied by changing the exponent a in the age- 
and connectivity- dependent attachment probability P(k, a) oc (1 + k)e~ aa . This model 
has a discontinuous transition in the relative size of the largest connected component: 
below a critical value a c there is an extensive number of topologically tree-like clusters of 
connected nodes and at a first-order critical point a giant cluster of linearly connected 
nodes emerges. We analyzed the properties of this phase transition with numerical 
simulations, finding that fluctuations of the order parameter scale linearly with system 
size and that, close to the critical point, the inverse of a characteristic length scale 
vanishes linearly with the distance from the critical point, suggesting that the transition 
is in the universality class of Id percolation. 
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